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ABSTRACT 


A  summary  is  given  of  basic  formulas  for  the  guiding  center 
motion  of  energetic  charged  particles  trapped  in  a  dipolar  mag¬ 
netic  field.  These  formulas  for  longitudinal  drift  rates, 
latitudinal  bounce  periods,  equatorial  gyroradii,  and  equatorial 
gyroperiods  sore  then  stated  in  convenient  numerical  form  for 
electrons  and  protons  as  functions  of  kinetic  energy  E,  magnetic^  q 

shell  parameter  L,  and  equatorial  pitch  angle/a^  for  a  slightly 
simplified  model  of  the  observed  magnetic  field  of  Saturn.  To 
aid  in  the  study  of  the  interaction  of  charged  particles  with  the 
rings  and  inner  satellites  of  Saturn,  additional  formulas  are 
given  for  the  time  interval  betveen  successive  encounters  of 
charged  particles  with  a  satellite  in  a  circular  prograde  orbit 
and  for  the  energies  of  electrons  whose  longitudinal  angular 
velocity  is  resonant,  or  synchronous,  with  the . Keplerian 
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INTRODUCTION 


The  purpose  of  this  short  paper  is  to  summarize  useful 

parameters  which  characterize  the  motion  of  charged  particles  in 

Saturn's  inner  magnetosphere.  Such  motion  is  a  combination  of 

gyration  about  a  magnetic  field  line,  latitudinal  bounce  along 

the  field  line,  and  longitudinal  drift  perpendicular  to  the 

field  line.  Hie  longitudinal  drift  arises  from  the  presence  of 

a  (corotational)  electric  field  and  from  the  curvature  and 

gradient  of  the  magnetic  field.  The  adopted  magnetic  field 
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model  is  that  of  a  point  dipole  of  moment  M  »  0.20  gauss  R 

s 

OD  X 

=  4.32  x  10  °  gauss  cat  (equatorial  radius  of  the  planet 
=  1  R  -  6  x  10^  cm),  located  at  the  geometrical  center  of  the 
planet  with  M  parallel  to  and  in  the  same  sense  as  Q,  the 
rotational  angular  velocity  of  the  planet.  This  model  is  a 
slight  simplification  of  the  observed  field  [Smith  et  al. , 

1980 ;  Acuna  and  Ness,  1980].  Parameters  related  to  the 
motion  of  electrons  and  protons  are  presented  as  functions  of 
kinetic  energy  E,  magnetic  shell  parameter  L,  and  equatorial 
pitch  angle  aQ.  We  adopt  the  convention  that  the  longitudinal 
angular  velocity  is  positive  for  eastward  drift  (same  sense  as 
the  rotation  of  the  planet)  and  negative  for  westward  drift. 


Unless  otherwise  stated,  all  basic  formulas  we  in  gaussian  c.g.s. 
units  and  all  numerical  formulas  give  angular  velocities  in  radian 
s  ^  distances  in  cm  and  periods  in  seconds  when  kinetic  energy  E 
is  expressed  in  MeV.  A  short  table  of  exemplary  values  is 
included.  Extensive  tables  are  contained  in  University  of 
Iowa  80-25  of  July  1980,  available  on  request  to  the  authors. 


SYMBOLS  AMD  CONSTANTS 

10  -1 

c  =  speed  of  light  *  2.9S7925  x  10  cm  s 

me2  *  rest  energy  of  the  particle 
*  0.511  MeV  for  electron. 

«  938.3  MeV  for  proton 

E  »  kinetic  energy  of  particle 

v  «  speed  of  particle 

P  »  v/c 

7  a  (1  -  P2)"^ 

q  =  charge  of  particle  *  4.80325  X  10  e.s.u. 
per  unit  charge 

R  ■  equatorial  radius  of  planet;  1  Rg  *  6  x  10  cm 
for  Saturn 

B  ■  equatorial  surface  field  of  planetary  dipole  moment 
o 

a  o.20  gauss  for  Saturn 
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M  *  dipole  moment  of  planet  »  4. 52  x  10  gauss  cm 
for  Saturn 

r  *  radial  distance  from  center  of  the  planet 

ao  *  equatorial  pitch,  angle  of  a  particle 

L  =  equatorial  crossing  distance  of  a  dipolar 
magnetic  field  line  in  units  of  R 

\  a  magnetic  latitude 

\  »  magnetic  latitude  of  a  particle's  mirror  point 

m 

B(\),  B^) 

■  magnetic  field  strength  at  \f  \a 
GM  *  gravitational  constant  of  planet 

pp  x  _p 

a  3.79311  x  10  cm  s  for  Saturn 

J_  a  second  zonal  harmonic  coefficient  of  planet's 

gravitational  potential  =>  1667  X  10~'  for  Saturn 

a  a  semimajor  axis  of  a  satellite  orbit 


Other  symbols  are  defined  as  they  appear  in  the  text. 
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GRADIENT  AND  CURVATURE  DRIFT 


The  bounce-averaged  longitudinal  drift  of  the  guiding  center 
of  a  particle  trapped  in  a  dipolar  magnetic  field  due  to  the 
curvature  and  gradient  of  the  field  has  been  studied  by  Hamlin 
et  al.  [1961],  Lew  [1961]  and  subsequent  authors.  Following 
Hamlin  et  al.  [1961],  Lew  [1961]  found  that  the  bounce -averaged 
angular  velocity  can  be  expressed  as 


“D 


SmcVyL 

a*/ 


where 


(1) 


(1  -  sin1*  \)  . 

(1+3  sin2  \)^ 

(1  +  3  alng  \)* 
Cl  -  BOO/B^)]^ 


2-B(A)/B(Xm) 

[1-B(\)/B(\a)]i 

cos\d\  . 


cos  \  d  \ 


The  integrals  are  carried  out  over  a  complete  latitudinal  bounce 
period.  For  Saturn,  Equation  (l)  becomes 


■»  *  *  ^ x  10'5  LE  (^2?)  ( 5)  • 


In  (2)  the  +  sign  applies  to  protons  and  the  -  sign  to  electrons. 

Lew  [196l]  further  presented  an  analytical  approximation  to  the 

ratio  F(\  )/G(X  )  which  is  good  to  within  one  part  in  10^. 
m  m 


(f/G)-1  »  [A+B  sin2  X  -  C  exp  (-k  sin2  X  )] 

'  m  m 


1.04675 

0.45555 

0.04675 

6.54568 


Values  of  X  and  F/G  for  several  values  of  the  particle's 
m  ' 

equatorial  pitch  angle  ctQ  are  given  in  Table  1. 

It  is  essential  to  note  that  is  the  longitudinal 
angular  velocity  of  particle  drift  in  the  reference  frame  rotat¬ 
ing  with  the  planet. 


COROTATIONAL  DRIFT 


In  addition  to  longitudinal  drift  due  to  the  curvature  and 
gradient  of  the  magnetic  field,  charged  particles  also  drift  in 
the  presence  of  an  electric  field  perpendicular  to  the  magnetic 
field  [Jackson,  1962].  Referenced  to  an  inertial  coordinate  system, 
an  electric  field  of  this  nature  exists  in  the  magnetosphere  of  a 
rotating  planet  if  the  planetary  ionosphere  is  electrically  coupled 
to  the  magnetosphere,  as  is  the  ease  in  the  presence  of  the  iono¬ 
spheric  and  magnetospheric  plasmas  which  make  the  electrical 
conductivity  parallel  to  magnetic  field  lines  essentially  infinite. 
Davis  [19*46]  showed  that  charge  flows  in  such  a  plasma  surrounding 
a  rotating,  conducting  magnetized  sphere  until  the  electric  field 

E  =  -i(nxr)xB  (4) 

is  established  in  the  plasma.  This  electric  field  results  in  a 
longitudinal  drift  angular  velocity  equal  to  the  planet's  angular 
velocity  fi  and  in  the  same  sense  [Birmingham  and  Northrop,  1979]* 

The  value  adopted  herein  for  0,  the  angular  velocity  of  Saturn's 
rotation,  is  1.637  X  10 radian  s”1,  which  corresponds  to  the 
rotational  period  10  hours  39*9  ±  0.3  minutes  determined  recently 
from  radio  observations  of  Saturn  by  Voyager  instruments  [Kaiser 
et  al.,  I960]. 
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ANGUIAR  VELOCITY  IN  AN  INERTIAL  FRAME 


The  corotational  angular  velocity  is  additive  to  the 
magnetic  gradient  and  curvature  drift  angular  velocity.  The 
algebraic  sum  of  these  t*,o  angular  velocities  is  the  longitudinal 
angular  velocity  of  a  particle  relative  to  an  inertial  frame. 

Thus 


<i)j  =  H  +  ujp  . 


(5) 


KEPLER  IAN  ORBITAL  ANGUIAR  VELOCITY 


From  Aksnes  [1$773>  the  mean  angular  velocity  of  a  neutral 
body  in  gravitational  orbit  about  Saturn  is 


GM 

t  [ 

3 

\a  i 

3J2  R 

“ST 


(6) 


With  a  in  units  of  R, 


%  (i.isSLSJ wf  ] . 


(7) 
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For  r  =  1.871  Rg,  0  =  u>k;  for  r  >  1.871  R0,  0  >  u^;  and 
for  r  <  1.871  Rg>  fi  <  id^  . 

TIME  INTERVAL  BETWEEN  SUCCESSIVE  ENCOUNTERS 
OF  A  CHARGED  BARTICLE  WITH  A  SATELLITE 

The  angular  velocity  of  a  charged  particle  relative  to  a 
satellite  in  a  circular  prograde  orbit  in  the  equatorial  plane 
at  the  same  radial  distance  is 

(»i  -  \)  =  (n  +  ^  -  (8) 

and  the  interval  between  encounters  is 

TE  =  2tt/(K  "  «kl)  •  (9) 

Note  that  (u)^  -  ux^)  may  be  positive  (eastward),  zero  (resonant), 
or  negative  (westward). 

RESONANT  ENERGY 

The  value  of  E  such  that  T^  is  infinite  is  called  the 
resonant  (or  synchronous)  energy  Ep  and  is  given  by 
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fi  +  (Op  -  =  0, 


(10) 


a  quadratic  equation  in  Ep,  whose  solution  for  electrons  is 


Ep  =  0.5  [E'  -  1.022  +  (E'2  +  1.0445)^]  , 


(11) 


where 


E'  =  4.800  x  101*  (fi  -  m^)  /  (LF/G)  .  (12) 

Only  for  r  <  1.871  R_  can  protons  be  resonant  with  a  satellite; 

s 

inasmuch  as  this  dividing  radius  is  near  the  outer  edge  of  Ring  B 
there  is  no  evident  interest  in  proton  resonances  in  Saturn's 
magneto  sphere . 


BOUNCE  PERIOD 


The  latitudinal  bounce  period  of  a  particle  oscillating 
between  mirror  points  on  a  given  magnetic  field  line  is 

%  ■  f  ^  -  <«> 

where  the  integral  is  carried  out  over  arc  length  s  along  the 
field  line  for  a  complete  (round  trip)  bounce  between  mirror 
points.  For  a  dipole  field, 
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a) 


g 


y  me 


qBe  . 
2 

E  tie 


(19) 


In  the  equatorial  plane  of  the  Saturnian  dipolar  magnetic  field. 


S0 


g 


1.798  X 10° 
L3  (E  +  me2) 


(20) 


The  corresponding  equatorial  gyroperiod  T  *  2tt/u>  . 

6  6 


T  -  3.495  x  10"6  L3  (E  +  me2) 

s 


(21) 


EQUATORIAL  GYRORADIUS 


The  equatorial  gyroradius  (radius  of  the  cylindrical 
surface  on  which  the  particle's  helical  trajectory  lies). 


r 


g 


S  c  sin  a 

o 


u> 


g 


(22) 


or 


1.667  X  10^  L3  sin  a 


[E(E  +  2mc2))^ 


(23) 


APPLICABILITY  OF  FORMULAS 


The  formulas  given  in  the  foregoing  sections  are  useful 
for  any  values  of  E,  L,  and  0^  for  which  the  assumptions  are 
applicable.  For  L  ^  7*  the  adopted  model  of  Saturn's  magnetic 
field  becomes  progressively  less  accurate  because  of  magneto- 
spheric  current  systems  and  for  L  i  13  the  formulas  are  of  very 
limited  value.  Also,  they  are  not  applicable  to  particles  whose 
energies  exceed  the  stable  trapping  limit  [Van  Allen,  1962]. 


L  VALUES  OF  SPECIAL  INTEREST 

A  particular  purpose  of  this  paper  is  to  aid  the  study  of 
the  interaction  of  trapped  particles  with  the  rings  and  inner 
satellites  of  Saturn.  L  values  of  special  interest  are  given  in 
Table  2.  In  refined  considerations,  a  full  set  of  orbital 
elements,  including  eccentricity  and  inclination,  is  required 
[Explanatory  Supplement,  1961]  as  are  diameters  of  the  satellites 
[Cruikshank,  1978]. 
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EXEMPIARY  TABLE 

In  Table  3,  we  give  exemplary  values  of  the  parameters 
discussed  herein  for  L  =  3.092.  Note  that  in  this  table  we 
depart  from  the  system  of  units  in  which  the  numerical  formulas 
are  given. 
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Table  1 


a 

0 

A. 

m 

F/G 

H 

90° 

0?0 

1.000 

0.74o 

80 

4.7 

0.995 

0.747 

70 

9-6 

0.980 

0.769 

6o 

14.7 

0.957 

0.805 

50 

20.2 

0.927 

0.855 

4o 

26.3 

O.891 

0.918 

30 

33.2 

0.851 

0.994 

20 

4l.4 

0.805 

1.083 

10 

52.5 

0.751 

1.191 
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Table  2 


L  Values  of  Stjecial  Interest 

Satellite 

L 

References 

1979  S5/1979  S6 

2.347 

Van  Allen  et  al.  [1960b] 

1979  S2/1979  s4 

2.528 

Van  Allen  et  al.  [1960b] 

1979  S3 

2.840 

Van  Allen  et  al.  [1960a;  1960b] 

SI  Mimas 

3.092 

Null  and  Lieske  [i960]; 
Explanatory  Supplement  to  the 
Astronomical  Ephemeris  and 
Nautical  Almanac  [1961] 

SII  Enceladua 

3.968 

Explanatory  Supplement  . . . 

[1961] 

Sill  Tethys 

4.913 

Explanatory  Supplement  . . . 

[1961] 

8IV  Dione 

6.292 

Explanatory  Supplement  ... 

[1961] 

Explanatory  Supplement  . . . 

[1961] 


8.787 


Values  of  Trapped  Particle  Parameters  for 
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